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ABSTRACT 

We outline our methods for obtaining high precision mass profiles, combining independent weak-lensing dis- 
tortion, magnification, and strong-lensing measurements. For massive clusters the strong and weak lensing 
regimes contribute equal logarithmic coverage of the radial profile. The utility of high-quality data is limited 
by the cosmic noise from large scale structure along the line of sight. This noise is overcome when stacking 
clusters, as too are the effects of cluster asphericity and substructure, permitting a stringent test of theoreti- 
cal models. We derive a mean radial mass profile of four similar mass clusters of high-quality Hubble Space 
Telescope and Subaru images, in the range R = 40kpc/i"' to 2800kpc/i"', where the inner radial boundary is 
sufficiently large to avoid smoothing from miscentering effects. The stacked mass profile is detected at 58(t 
significance over the entire radial range, with the contribution from the cosmic noise included. We show that 
the projected mass profile has a continuously steepening gradient out to beyond the virial radius, in remarkably 
good agreement with the standard Navarro-Frenk- White form predicted for the family of CDM-dominated ha- 
los in gravitational equilibrium. The central slope is constrained to lie in the range, -dlnp/dlnr = 0.89;|;o 39. 
The mean concentration is Cvir = 7.68;!;0 4q (^tMvir = L54;|;q x lO'^Mo/i"'), which is high for relaxed, high- 
mass clusters, but consistent with ACDM when a sizable projection bias estimated from A^-body simulations is 
considered. This possible tension will be more definitively explored with new cluster surveys, such as CLASH, 
LoCuSS, Subaru HSC, and XXM-XXL, to construct the Cvii-Mvir relation over a wider mass range. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: general — gravitational 
lensing: weak — gravitational lensing: strong 



L INTRODUCTION 

Clusters of galaxies represent the largest gravitationally- 
bound objects in the universe, which contain a wealth of 
astrophysical and cosmological information, related to the 
nature of dark matter, primordial density perturbations, and 
the emergence of structure over cosmic time. Observational 
constraints on the properties and evolution of clusters pro- 
vide independent and fundamental tests of any viable cos- 
mology, structure formation scenario, and possible modifica- 
tions of the laws of gravity, complementing large-scale cos- 
mic microwave background and galaxy clustering measure- 
ments (e.g., Komatsu et al. 201 1; Percival et al. 2010). 

A key ingredient of cluster-based cosmological tests is the 
mass and internal mass distribution of clusters. In this con- 
text, the current cosmological paradigm of structure forma- 
tion, the standard A cold (i.e., non relativistic) dark matter 
(ACDM, hereafter) model, provides observationally testable 
predictions for CDM-dominated halos over a large dynamical 
range in density and radius. Unlike galaxies where substantial 
baryonic cooling is present, massive clusters are not expected 
to be significantly affected by gas cooling (e.g., Blumenthal 
et al. 1986; Broadhurst & Barkana 2008). This is because 
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the majority of baryons (^^ 80%) in massive clusters comprise 
a hot. X-ray emitting diffuse intracluster medium (hereafter 
ICM), in which the high temperature and low density prevent 
efficient cooling and gas contraction, and hence the gas pres- 
sure roughly traces the gravitational potential produced by the 
dominant dark matter (see Kawaharada et al. 2010; Molnar 
et al. 2010). The ICM represents only a minor fraction of 
the total mass near the centers of clusters (Lemze et al. 2008 ; 
Umetsu et al. 2009). Consequently, for clusters in a state of 
quasi equilibrium, the form of their total mass profiles reflects 
closely the distribution of dark matter (Mead et al. 2010). 

High-resolution A^-body simulations of collisionless CDM 
exhibit an approximately "universal" form for the spherically- 
averaged density profile of virialized dark matter halos 
(Navarro et al. 1997, NEW hereafter), with some intrin- 
sic variance in the mass assembly histories of individual 
halos (Jing & Suto 2000; Tasitsiomi et al. 2004; Navarro 
et al. 2010). The predicted logarithmic gradient 73d('") = 
-dlnp/dlnr of the NEW form flattens progressively toward 
the center of mass, with a central cusp slope flatter than 
a purely isothermal structure (730 = 2) interior to the in- 
ner characteristic radius rs{ ^ 300 kpc/;"' for cluster-sized ha- 
los), providing a distinctive prediction for the empirical form 
of CDM halos in gravitational equilibrium. A useful index 
of the degree of concentration is Cvir = rya/rs, which com- 
pares the virial radius r^h- to the characteristic radius rj of the 
Navarro-Frenk- White (NFW, hereafter) profile. This empiri- 
cal NEW profile is characterized by the total mass within the 
virial radius, Mvir, and the halo concentration Cvir- Theoretical 
progress has been made in understanding of the form of this 
profile in terms of the dynamical structure using both numer- 
ical and analytical approaches (Taylor & Navarro 2001; Lapi 
& Cavaliere 2009; Navarro et al. 2010), though we must cur- 
rently rely on the quality of A^-body simulations when making 
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comparisons with CDM-based predictions for cluster mass 
profiles. 

In the context of standard hierarchical clustering models, 
the halo concentration should decline with increasing halo 
mass as dark matter halos that are more massive collapse later 
when the mean background density of the universe is corre- 
spondingly lower (Bullock et al. 2001; Zhao et al. 2003; Neto 
et al. 2007). This prediction for the halo mass-concentration 
relation and its evolution has been established thoroughly 
with detailed simulations (e.g., Navarro et al. 1997; Bullock 
et al. 2001; Neto et al. 2007; Duffy et al. 2008; Klypin et al. 

2010) , although sizable scatter around the mean relation is 
present due partly to variations in formation epoch of halos 
(Wechsler et al. 2002; Neto et al. 2007; Zhao et al. 2009). 
Massive clusters are of particular interest in this context be- 
cause they are predicted to have a relatively shallow mass pro- 
file with a pronounced radial curvature. Gravitational lensing 
of background galaxies offers a robust way of directly obtain- 
ing the mass distribution of galaxy clusters (see Bartelmann & 
Schneider 2001; Umetsu 2010, and references therein) with- 
out requiring any assumptions on the dynamical and physi- 
cal state of the system (Clowe et al. 2006; Okabe & Umetsu 
2008). A detailed examination of this fundamental predic- 
tion has been the focus of our preceding work (Broadhurst 
et al. 2005b; Medezinski et al. 2007; Broadhurst et al. 2008; 
Umetsu & Broadhurst 2008; Umetsu et al. 2009; Lemze et al. 
2009; Umetsu et al. 2010, 201 1). 

Systematic cluster lensing surveys are in progress to ob- 
tain mass profiles of representative clusters over a wide range 
of radius by combining high-quaUty strong and weak lens- 
ing data. Deep multicolor images of massive cluster cores 
from Advanced Camera for Surveys (ACS) observations with 
the Hubble Space Telescope (HST) allow us to identify many 
sets of multiple images spanning a wide range of redshifts 
for detailed strong-lens modeling (e.g., Broadhurst et al. 
2005a; Zitrin et al. 2009, 2010, 2011a,c,d). The wide-field 
prime-focus cameras of Subaru and CFHT have been rou- 
tinely producing data of sufficient quality to obtain accurate 
measurements of the weak-lensing signal, providing model- 
independent cluster mass profiles out to beyond the virial 
radius (e.g., Broadhurst et al. 2005b, 2008; Limousin et al. 
2007; Umetsu & Broadhurst 2008; Umetsu et al. 2009, 2010, 
2011; Coe et al. 2010). Our earlier work has demonstrated 
that without adequate color information, the weak-lensing 
signal can be heavily diluted particularly toward the cluster 
center by the presence of unlensed cluster members, lead- 
ing to biased cluster mass profile measurements with under- 
estimated concentrations and internal inconsistency, with the 
weak-lensing based profile underpredicting the observed Ein- 
stein radius (Broadhurst et al. 2005b; Umetsu & Broadhurst 
2008; Medezinski et al. 2010). 

Careful lensing work on individual clusters has shown that 
full mass profiles constructed from combined strong and weak 
lensing measurements show a continuous steepening radial 
trend consistent with the predicted form for the family of col- 
lisionless CDM halos (Gavazzi et al. 2003; Broadhurst et al. 
2005b, 2008; Umetsu & Broadhurst 2008; Umetsu et al. 2010, 

201 1) . Intriguingly these initial results from combined strong 
and weak lensing measurements reveal a relatively high de- 
gree of halo concentration in lensing clusters (e.g., Gavazzi 
et al. 2003; Kneib et al. 2003; Broadhurst et al. 2008; Oguri 
et al. 2009; Zitrin et al. 2011d), lying well above the mass- 
concentration relation for cluster-sized halos predicted by the 
ACDM model, despite careful attempts to correct for potential 



projection and selection biases inherent to lensing (Hennawi 
et al. 2007; Meneghetti et al. 2010b). This apparent overcon- 
centration of lensing clusters is also indicated by the gener- 
ally large Einstein radii determined from strong-lensing data 
(Broadhurst & Barkana 2008; Meneghetti et al. 2010a; Zitrin 
etal. 2011a). 

In this paper we explore in greater depth the utility of 
high-quality lensing data for obtaining highest-precision clus- 
ter mass profiles by combining all possible lensing informa- 
tion available in the cluster regime. This extends our recent 
weak-lensing work by Umetsu et al. (201 1), where a Bayesian 
method was developed for a direct reconstruction of the pro- 
jected cluster mass profile from complementary weak-lensing 
distortion and magnification effects (Umetsu & Broadhurst 
2008), the combination of which can be used to unambigu- 
ously determine the absolute mass normalization. For a mas- 
sive cluster acting as a super-critical lens, the strong and weak 
lensing regimes contribute equal logarithmic coverage of the 
radial profile (Umetsu et al. 20 11), so that here we concen- 
trate on those clusters for which we have high-quahty data 
in both these regimes. The high quality of our data is such 
that we have now become significantly limited by the cos- 
mic noise from large scale structure behind the cluster center, 
where magnified sources lie at greater distances. This noise 
is correlated between radial bins, and so can be overcome by 
stacking clusters, along independent sight Unes. Stacking also 
helps average over the effects of cluster asphericity and sub- 
structure, (Mandelbaum et al. 2006; Johnston et al. 2007; Ok- 
abe et al. 2010; Umetsu et al. 201 1), allowing a tighter com- 
parison of the averaged profile with theoretical models. Our 
aim here is to develop and apply comprehensive methods to 
a sample of four similarly high-mass lensing clusters (A 168 9, 
A 1703, A370, and C10024H-17), for which we have previously 
identified multiply-lensed images and measured weak mag- 
nification and distortion effects from deep HST and Subaru 
observations (Broadhurst et al. 2005a; Umetsu & Broadhurst 
2008; Umetsu et al. 2010; Zitrin et al. 2010; Medezinski et al. 
2010, 201 1; Umetsu et al. 201 1). 

The paper is organized as follows. In Section 2 we briefly 
sunmiarize the basis of cluster gravitational lensing. In Sec- 
tion 3 we outline our comprehensive methods for obtaining 
projected cluster mass profiles from weak-lensing distortion, 
magnification, and strong-lensing measurements. In Section 
4 we apply our methodology to deep HST and Subaru obser- 
vations of four massive clusters to derive a mean radial mass 
profile for the entire clusters, demonstrating how stacking the 
weak and strong lensing signals improves upon the statistical 
precision of the mass profile determination; we then examine 
the radial dependence of the stacked cluster mass profile. Fi- 
nally, we discuss our results and conclusions in § 5. Through- 
out this paper, we adopt a concordance ACDM cosmology 
with n,„ = 0.3, = 0.7, and h = //o/(100kms-'Mpc-i) = 
0.7, unless otherwise noted. 

2. BASIS OF CLUSTER LENSING 

The gravitational deflection of light rays by a cluster can be 
described by the thin lens equation, which relates the angular 
position of a lensed image 9 to the angular position of the 
intrinsic source /3 as 

/3 = 0-VV, (1) 

where a = 'Vtp(O) is the deflection field, and tl^(9) 
is the effective lensing potential, which is defined by 
the two-dimensional Poisson equation as Atjj{6) = 2K{d) 
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with the lensing convergence k given as a source term. 
This equation can be readily inverted to yield: tp(6) = 

ijd^e' A-i (61, e'MO') = (i/tt) _/ d^e'in\e-e'\ K{e'), so that 

the deflection field is expressed in terms of k, as 
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For gravitational lensing in the cluster regime (e.g., Umetsu 
2010), K is expressed as k(0) = E^jS(0), namely the pro- 
jected mass density S(0) in units of tiie critical surface mass 
density for gravitational lensing, defined as 



47rGD 



■r'; ^(z.) = max 



DlsiZs) 



(3) 



where D^, Di, and D/j are the proper angular diameter dis- 
tances from the observer to the source, from the observer to 
the lens, and from the lens to the source, respectively, and 
P = Dis/Ds is the angular-diameter distance ratio associated 
with the population of background sources. 

The deformation of the image for a background source can 
be described by the Jacobian matrix Aa/B = {df3/d6)oc^ = 
5ai3-i'.a0 (a,/3 = 1,2) of the lens mapping, where Sap is 
Kronecker's delta. The real, symmetric Jacobian Aajj can 
be decomposed as Aap ' (1 -'«)<5a/J~r„^, where r„^(0) = 
{dadp — Safi^'^/2)'tjj{6) is the trace-free, syrmnetric shear ma- 
trix. 
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with 7q being the components of spin-2 complex gravita- 
tional shear 7 := 71 + 1^2- In the strict weak-lensing limit 
where k, I7I <C 1, T^p induces a quadrupole anisotropy of 
the background image, which can be observed from elliptici- 
ties of background galaxy images. Given an arbitrary circular 
loop of radius on the sky, the average tangential shear 7+(i?) 
around the loop satisfies the following identity (e.g.. Kaiser 
1995): 

7+(^) = K«^)-/^(i?), (5) 

where K(t?) is the azimuthal average of Kid) around the loop, 
and k{< t?) is the average convergence within the loop. 

The local area distortion due to gravitational lensing, or 
magnification, is given by the inverse Jacobian determinant. 
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which can influence the observed surface density of back- 
ground sources, expanding the area of sky, and enhancing the 
observed flux of background sources (Broadhurst et al. 1995). 
The former effect reduces the effective observing area in the 
source plane, decreasing the number of background sources 
per sohd angle; on the other hand, the latter effect ampli- 
fies the flux of background sources, increasing the number 
of sources above the limiting flux. The net effect is known as 
magnification bias and depends on the intrinsic slope of the 
luminosity function of background sources. 

* Throughout the paper we assume in our weak lensing analysis that the 
angular size of background galaxy images is sufficiently small compared to 
the scale over which the underlying lensing fields vary, so that the higher- 
order weak lensing effects, such as flexion, can be safely neglected; see, e.g., 
Goldberg & Bacon (2005); Okura et al. (2007, 2008). 



In general, the observable quantity for quadrupole weak 
lensing is not the gravitational shear 7 but the complex re- 
duced shear, 

7(0) 

1 - K{U) 

in the subcritical regime where det.4. > (or 1 /g* in the nega- 
tive parity region with det^A < 0). The spin-2 reduced shear g 
is invariant under the following global linear transformation: 



k(6/)^ Ak(6/)+1-A, 7(0)^A7(6») 



(8) 



with an arbitrary scalar constant X (Schneider & Seitz 
1995). This transformation is equivalent to scaling the Jaco- 
bian matrix Aid) with A, Aid) XAiO), and hence leaves the 
critical curves detAiO) = invariant. Furthermore, the curve 
niO) = 1, on which the gravitational distortions disappear, is 
left invariant under the transformation (8). 

This mass-sheet degeneracy can be unambiguously broken 
by measuring the magnification effects, because the magni- 
fication fj, transforms under the invariance transformation (8) 
as 



fxiO) ^ AV(0). 



(9) 



In practice, the lens magnification /i can be measured from 
characteristic variations in the number density of background 
galaxies due to magnification bias (Broadhurst et al. 1995; 
Umetsu etal. 2011) as 

n^(0) = noMW"^\ (10) 

where «o = dNQi< nicud/dfl is the unlensed number den- 
sity of background sources for a given magnitude cutoff 
Wcut, approximated locally as a power-law cut with slope 
s = dlogiQA'o(< m)/dm is > 0). In the strict weak-lensing 
limit, the magnification bias is dn^/no « (5s-2)k. For red 
background galaxies the intrinsic count slope s at faint mag- 
nitudes is relatively flat, ^ 0.1, so that a net count depletion 
results (Broadhurst et al. 2005b; Umetsu & Broadhurst 2008; 
Umetsu et al. 2010, 2011). On the other hand, the faint blue 
background population tends to have a steeper intrinsic count 
slope close to the lensing invariant slope is = 0.4). Alterna- 
tively, the constant A can be determined such that the mean 
K averaged over the outermost cluster region vanishes, if a 
sufficiently wide sky coverage is available.^ 

3. CLUSTER LENSING METHODOLOGY 

In this section we outline our methods for obtaining cluster 
mass profiles in a continuous radial coverage from the cen- 
tral region to beyond the virial radius, by combining inde- 
pendent weak-lensing distortion, magnification, and strong- 
lensing measurements. 

3.1. Cluster Weak Lensing 

The relation between observable distortion ig) and underly- 
ing convergence in) is nonlocal, and the convergence derived 
from distortion data alone suffers from a mass-sheet degener- 
acy (§2). However, by combining the complementary distor- 
tion ig) and magnification (/i) measurements the convergence 
can be obtained unambiguously with the correct mass normal- 
ization. 

We construct a discrete convergence profile in the weak- 
lensing regime from observable lens distortion and magnifica- 
tion profiles, g+iO) = 7+(6i)/[l - k(6')] and n^(6») = no/J-ie)^-^'-^ 

' Or, one may constrain the constant A such that the enclosed mass within 
a certain aperture is consistent with cluster mass estimates from some other 
observations (e.g., Umetsu & Futamase 2000). 
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(see Section 3 and Appendix B of Umetsu et al. 2011, for 
details of weak-lensing profile measurements), following the 
Bayesian prescription given by Umetsu et al. (2011). The 
Bayesian approach allows for a full parameter-space extrac- 
tion of model and caUbration parameters. A proper Bayesian 
statistical analysis is of particular importance to explore the 
entire parameter space and investigate the parameter degen- 
eracies, arising in part from the mass-sheet degeneracy. 

In the Bayesian framework, we sample from the posterior 
probabiUty density function (PDF) of the underlying signal s 
given the data d, P(s\d). Expectation values of any statistic of 
the signal s shall converge to the expectation values of the a 
posteriori marginalized PDF, P(s\d). The variance co variance 
matrix C of s is obtained from the resulting posterior sample. 
In our problem, the signal s is a vector containing the dis- 
crete convergence profile, k, = K(9i) with / = 1,2, ..jA'"' in the 
weak-lensing regime (Oj > ^Ein), and the average convergence 
within the inner radial boundary of the weak lensing data, 

Kmin = k(< ^^n)' ^h^^ * = {^min, Kij^p ^^^^E Specified by 
(A'^' + 1) parameters. The Bayes' theorem states that 

P(s\d)(xP(s)P(d\s), (11) 

where £(s) = Pid\s) is the likelihood of the data given the 
model (s), and P(s) is the prior probability distribution for 
the model parameters. The C(s) function for combined weak 
lensing observations is given as a product of the two sepa- 
rate likelihoods, £wi = ^C^-C^, where Cg and £^ are the likeli- 
hood functions for distortion and magnification, respectively, 
as given in Umetsu et al. (2011). The log-UkeUhood for 
combined weak-lensing distortion and magnification observa- 
tions, {g+,i}fJi and {n^jjf^i , is given as 

-21n£„i = ^ + 2^ / , (12) 

i=l i=l "f^.i 

where {g+^i,hfj,^i) are the theoretical predictions for the corre- 
sponding observations; the errors cr+,i forg+_, (/= 1,2, ...jA'*') 
due primarily to the variance of the intrinsic source ellip- 
ticity distribution can be conservatively estimated from the 
data using bootstrap techniques; the errors a^^i for n^^i (i = 
1,2, ...,N^^) include both contributions from Poisson errors in 
the counts and contamination due to intrinsic clustering of red 
background galaxies (Umetsu et al. 201 1). 

For each parameter of the model s, we consider a simple flat 
prior with a lower bound of s = 0, that is, Kmin > and k, > 0. 
Additionally, we account for the calibration uncertainty in the 
observational parameters, such as the normalization and slope 
parameters (no,s) of the background counts and the relative 
lensing depth due to population-to-population variations be- 
tween the background samples used for the magnification and 
distortion measurements (see Umetsu et al. 2011). 

3.2. Cluster Strong Lensing 

We apply our well-tested approach to strong-lens modehng, 
which has previously uncovered large numbers of multiply- 
lensed galaxies in ACS images of many clusters, such as 
A1689 at z = 0.183 (Broadhurst et al. 2005a), C10024H-17 at 
z = 0.395 (Zitrin et al. 2009), 12 high-z MACS clusters (Zitrin 
et al. 201 la), MS 1358H-62 at z = 0.33 (Zitrin et al. 201 Ic), and 
A383 at z = 0.188 (Zitrin et al. 201 Id). Briefly, tiie basic as- 
sumption adopted is that mass approximately traces light, so 
that the photometry of the red cluster member galaxies is used 
as the starting point for our model. Cluster member galaxies 



are identified as lying close to the cluster sequence by HST 
multiband photometry. 

In the strong-lensing regime we approximate the large scale 
distribution of cluster mass by assigning a power-law mass 
profile to each cluster galaxy, the sum of which is then 
smoothed. The degree of smoothing (S) and the index of 
the power-law (q) are the most fundamental parameters deter- 
mining the cluster mass profile dominated by dark matter. A 
worthwhile improvement in fitting the location of the lensed 
images is generally found by expanding to first order the grav- 
itational potential of this smooth component, equivalent to a 
coherent extemal shear F^"^ (a,/3 = 1,2) describing the over- 
all matter eUipticity. The direction (/)ex of the spin-2 external 
shear F^^ and its amplitude |7ex| are free parameters, allow- 
ing for some flexibility in the relation between the distribution 
of dark matter and the distribution of galaxies, which cannot 
be expected to trace each other in detail. 

The total deflection field aiO) = J^jajiO) = 

(I]-4/7r)/c/26)'(6/ - e')/\e - 6'\^Y.j^0) consists of 
the galaxy component agai(0), scaled by a factor K, the 
smooth cluster dark-matter component aDM(^), scaled by 
{\-K), and the extemal-shear component aex(^), 

aid) = KcXgdO) + (1 -K)cxDMie) + cx^iO), (13) 

where Q!ex,a(^) = (^'"^)ap^^i^ with AO being the displace- 
ment vector of the angular position 6 with respect to a fiducial 
reference position. The overall normalization Af of the model 
and the relative scaling K of the smooth dark matter compo- 
nent versus the galaxy contribution bring the total number of 
free parameters in the model to 6. This approach to strong 
lensing is sufficient to accurately predict the locations and in- 
ternal structure of multiple images, since in practice the num- 
ber of multiple images uncovered readily exceeds the number 
of free parameters, so that the fit is fully constrained. 

We use this preliminary model to delens the more obvious 
lensed galaxies back to the source plane by subtracting the de- 
rived deflection field. We then relens the source plane in order 
to predict the detailed appearance and location of additional 
counter images, which may then be identified in the data by 
morphology, internal structure and color. The best-fit strong- 
lensing model is assessed by minimizing the value in the 
image plane: 

2 ^ [0,-0,<g,5,AA,g,r")]2 

Xsi = 2^ ^2 ' (14) 

i ' 

where i runs over all lensed images, 6i(q,S,J\f,T^^) is the 
position given by the model, 0, is the observed image posi- 
tion, and (J, is the positional measurement error. For each 
model parameter, we estimate the la uncertainty by Ax^ = 
X^~xLin = 1 in the six-parameter space. The uncertain- 
ties for the E(0) field and the T,{9) profile are estimated by 
propagating the errors on the strong-lens model parameters, 

3.3. Combining Weak and Strong Lensing 

We derive a full-radial mass profile on an individual clus- 
ter basis by combining independent weak and strong lensing 
data, which can be compared for consistency in the region 
of overlap. In order to obtain meaningful radial profiles, one 
must carefully define the center of the cluster It is often as- 
sumed that the cluster mass centroid coincides with the posi- 
tion of the brightest cluster galaxy (BCG), whereas the BCGs 
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can be offset from the mass centroids of the corresponding 
dark matter halos (Johnston et al. 2007; Oguri et al. 2010; 
Oguri & Takada 201 1). Umetsu et al. (201 1) adopted the loca- 
tion of the BCG as the cluster center in their one-dimensional 
profile analysis of five massive clusters. A small offset of typ- 
ically ^ 20 kpc = doff is found by Umetsu et al. (20 11) be- 
tween the BCG and the dark matter center of mass recovered 
from strong-lens modeUng (Section 3.2). In the following we 
will adopt the BCG position as the cluster center, and limit 
our analysis to radii greater than /?nun = 2Joff = 40 kpc h~^, be- 
yond which the cluster miscentering effects on the E profile 
are negligible (see Section 10 of Johnston et al. 2007). 

Having defined the cluster center, we can construct a joint 
discrete mass profile S = {E(/?,)}^j as a function of the pro- 
jected radius R = Dfi by combining the weak and strong lens- 
ing K profiles: S(7?;) = w~[^k{9i) (i = l,2,...,N), where w,- is 
the lensing efficiency function, or the inverse critical surface 
mass density, w, = (Scrit.;)"' = (47rG/c^)Z)//3,', Note, the ; de- 
pendence arises because strong and weak lensing profiles with 
different depths are combined together. To simplify the anal- 
ysis, we exclude the strong-lensing data points in the region 
of overlap (typically, 6'Ein i^O ^ 2^^Ein) as well as the central 
weak-lensing bin TCmin, when defining the joint S profile. 

The formulation thus far allows us to derive covariance 
matrices Clf^ of statistical measurement errors for individ- 
ual cluster K profiles. Here we take into account the effect 
of uncorrelated large scale structure projected along the line 
of sight on the error covariance matrix Cfj" as C = C^'^' + C'^^ 

where C^^ is given as (Schneider et al. 1998; Hoekstra 2003; 
Dalai et al. 2005; Hoekstra et al. 201 1 ; Oguri & Takada 201 1) 

c)f c'-iDjmmiej). (15) 

Here C'^'^(l) is the weak-lensing power spectrum as a function 
of angular multipole / evaluated for a given source population 
and a cosmology, and Jo(l9i) is the Bessel function of the first 
kind and order zero averaged over the ith annulus between 1 

and 0,.2(> ^i,i), given as 

= no M^na ^2 H2Jm,2)-lOi,iJm,i)] ■ (16) 

We will assume the concordance ACDM cosmological model 
of Komatsu et al. (2011) and use the fitting formula of Pea- 
cock & Dodds (1996) to compute the nonlinear mass power 
spectrum that enters in equation (15). 

3.4. Stacked Lensing Analysis 

The utility of high-quality data is ultimately limited by the 
cosmic noise from large scale structure along the line of sight, 
producing covariance between radial bins, particularly behind 
the cluster center, where magnified sources lie at greater dis- 
tances. This noise is correlated between radial bins, but can 
be overcome by stacking an ensemble of clusters along inde- 
pendent lines of sight. Stacking also helps average over the 
effects of cluster asphericity and substructure inherent in pro- 
jected lensing measurements. The statistical precision can be 
greatly improved by stacking together a number of clusters, 
especially on small angular scales (see Okabe et al. 2010), al- 
lowing a tighter comparison of the averaged profile with the- 
oretical models. 

With the full mass profiles of individual clusters from 
combined weak and strong lensing (Section 3.3), we can 
stack the clusters to produce an averaged radial mass pro- 
file. Here we re-evaluate the mass profiles of the individual 



clusters in M logarithmically-spaced radial bins in the range 
R = [^mini^maxl followiug the prescription given in Umetsu 
et al. (2011). Since the noise in different clusters is uncorre- 
lated, the mass profiles of individual clusters can be co-added 
according to (Umetsu et al. 201 1) 

{^)=(^w}j (^^W„S„j, (17) 

where the index n runs over all clusters, S„ is a vector con- 
taining the discrete surface mass density profile for the nth 
cluster, and yV„ is the window matrix defined as 

(Wn)ij^{C;')..{wMw„)j (18) 

with (Cn)ij and (w„), (i = 1,2, ...,M) being the full covariance 
matrix and the lensing efficiency function for the nth cluster, 
respectively. The error covariance matrix for the stacked mass 
profile (S) is obtained as 

4. APPLICATIONS: HUBBLE AND SUBARU OBSERVATIONS OF 
FOUR HIGH-MASS CLUSTERS 

4.1. Cluster Sample and Lensing Data 

Following the methodology outlined in Section 3, we an- 
alyze our consistent weak and strong lensing measurements 
presented in Umetsu et al. (2011) to examine the underlying 
projected mass profile S(/?) of a sample of four well-studied 
high-mass clusters (M ^ lO'^M©) at intermediate redshifts, 
A1689 (z = 0.183), A1703 (z = 0.281), A370 (z = 0.375), 
and C10024H-17 (z = 0.395)^". The massive clusters we have 
analyzed are well-known strong lensing clusters, displaying 
prominent strong-lensing features and large Einstein radii of 
^Ein ^ 30" (e.g., for a fiducial source redshift Zs ~ 2; Broad- 
hurst & Barkana 2008; Oguri & Blandford 2009; Zitrin et al. 
2011b). Table 1 gives a summary of the basic properties of 
the clusters in our sample. 

For these clusters, the central mass distributions 
(R ^ 200kpc h~^) have been recovered in detail by our 
strong-lensing analysis (Broadhurst et al. 2005a; Zitrin et al. 
2009, 2010; Umetsu et al. 2011) based on many sets of 
multiply-lensed images identified previously in very deep 
multicolor imaging with HST/ACS (e.g., Broadhurst et al. 
2005a; Limousin et al. 2008; Richard et al. 2009, 2010; 
Zitrin et al. 2009, 2010). Umetsu et al. (2011) developed 
and applied a Bayesian method to derive model-independent 
projected mass profiles for five high-mass clusters (including 
RXJ1347-11 in addition to the four clusters) from Subaru 
weak-lensing distortion and magnification measurements, the 
combination of which can unambiguously break the mass- 
sheet degeneracy inherent in any mass inversion method 
based solely on shape distortion data. It was shown that 
for the four clusters of the present sample our independent 
strong and weak lensing mass profiles are in full agreement 
in the region of overlap (R ~ 150kpc/i"'), and together can 
be well described by, within the noise, a generalized form 
of the NFW profile for CDM-dominated equilibrium halos. 
This motivates us to reexamine in detail the form of the radial 
mass profile for the entire clusters. 

Careful examination of lensing, X-ray, and dynamical data strongly sug- 
gest that C10024+17 is the results of a high-speed, line-of-sight collision of 
two massive clusters viewed approximately 2-3 Gyr after impact when the 
gravitational potential has had time to relax in the center, but before the gas 
has recovered (see Umetsu et al. 2010, and references therein) 
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4.2. Results 

Our weak and strong lensing data together cover a wide 
range of radius ranging typically from R ^ lOkpc/z"' to 2000- 
3500kpc/!"' (Umetsu et al. 2011), depending on the cluster 
redshift as limited by the field of view of Subaru/Suprime- 
Cam (34' x 27'). Table 1 hsts for each cluster the radial 
ranges R = [/?^„,/?lJ and of strong and weak 

lensing measurements, respectively, used to define a joint dis- 
crete mass profile S = {S(/?,)};Ip given in a total of A' radial 
bins spanning from T^^in = ^^in to ^max = Rmax- In Table 1, 
we also quote values of the total S/N ratio in our joint clus- 
ter mass profiles (S) obtained using the full covariance matrix 
C. We find that ignoring the cosmic noise contribution (equa- 
tion [15]) will underestimate the errors by ~ 30%^0%. To 
evaluate C'*** for strong-lensing observations, we projected the 
matter power spectrum out to a fiducial depth of Zs = 2, which 
is a typical source redshift of strongly-lensed arcs in clusters 
at intermediate redshifts. We used the estimated mean source 
redshifts given in Table 3 of Umetsu et al. (2011) for weak 
lensing. 

We show in the top panel of Figure 1 the resulting averaged 
radial mass profile (S(/?)) in M = 15 logarithmically-spaced 
bins with its statistical 1 cr uncertainty (given as the square root 
of the diagonal part of the full covariance matrix C), obtained 
by stacking the four clusters using equations (17) and (19). 
Note, no scaling has been applied to match the mass normal- 
izations between the four clusters, which span a relatively nar- 
row range in mass, 1 .3 < Mvir/(10'^Mo /;"') < 2.3 (see Table 
6 of Umetsu et al. 201 1). For our sample, we find a sensitivity- 
weighted average cluster redshift of (zi) ^ 0.32, which is 
fairly close to the simple average of z; = 0.31 due to the nar- 
row redshift coverage of our cluster sample. The stacked mass 
profile exhibits a smooth radial trend with a clear radial cur- 
vature over a wide range of radius from R = 40kpc/!~' to 
2800kpc/!"' w 1.4rvir, and is detected at a high significance 
level of 58(7, with the contribution from cosmic covariance 
included. Here the maximum radius for the stacking analysis 
represents approximately the average maximum radius (/?max) 
covered by our data. Also shown in Figure 1 is the cosmic 
noise contribution, which increases toward the cluster center. 
A noticeable increase of the stacked cosmic noise is seen at 
R ^ 150kpcft~\ within which the averaged profile is domi- 
nated by strong lensing measurements with greater depth. In 
the bottom panel of Figure 1 , we plot the logarithmic density 
slope 72d(^) ^ -dln{Tj) /dlnR of the stacked mass profile. 
The logarithmic gradient of the average profile shows a clear 
continuous steepening with increasing radius in projection. 

To quantify and characterize the averaged cluster mass dis- 
tribution, we compare the (S) profile with the physically and 
observationally motivated NFW model. Here we consider a 
generalized parametrization of the NFW model (gNFW, here- 
after) of the form (Zhao 1996; Jing & Suto 2000): 

P^''^ = r I ^on^I / ' (20) 
(r/rj)"(l + r/r,)J " 

where ps is the characteristic density, is the characteristic 
scale radius, and a represents the inner slope of the density 
profile. This reduces to the NFW model for a = 1. We in- 
troduce the radius r-2 at which the logarithmic slope of the 
density is isothermal, i.e., 730 = 2. For the gNFW profile, 
r_2 = (2-a)ri, and thus the corresponding concentration pa- 
rameter reduces to c_2 = r^-u/f-i = Cvii-/(2- a). We specify 
the gNFW model with the central cusp slope, a, the halo 
virial mass, Myir, and the concentration, c-2 = Cvir/(2-a). We 



employ the radial dependence of the gNFW lensing profiles 
given by Keeton (2001). 

First, when the central cusp slope is fixed to a = 1 (NFW), 
the best-fit model for the averaged (S) profile is obtained as 
Mvi,- = 1.54+!]:[i X lO'^Mo/i-i and c_2 = c^i, = 7.68^1];^^ with 
the minimized value (x^n) of 5.8 for 13 degrees of free- 
dom (dof), corresponding to a Q- value goodness-of-fit oiQ = 
0.952. This model yields an Einstein radius of = 39.9"^4 ^ 
for a fiducial source at z, = 3. The resulting best-fit NFW 
parameters from the stacked analysis are consistent with the 
respective sample weighted means of the individual NFW 
model fits obtained by Umetsu et al. (2011, Table 6): (Mvir) = 
1.44 ±0.11 X IO'^Mq/!-' and (cvir) = 7.76 ± 0.79. Next, 
when a is allowed to vary, a gNFW fit to (S) gives Myk = 
1.50!^:j^ X 1O>^M0 h-\ c-2 = and a = 0.89!;; ^ with 

xL/dof= 5.7/12 and Q = 0.931 (f?Em = 38.4"!i,^i at z, = 3), 
being consistent with a simple NFW model with a= 1. Thus 
the addition of the a parameter does not improve the fit 
substantially, as shown by the quoted and Q values (see 
also Zitrin et al. 201 Id). The two-dimensional marginalized 
constraints (68.3%, 95.4%, and 99.7% confidence levels) on 
(Mvir, ct) and (c_2, a) are shown in Figure 2. Finally, a force fit 
to the singular isothermal sphere (SIS) model (p cx r'^) yields 
a poor fit with x^n/dof = 78.5/14, so that the SIS model is 
strongly disfavored at 62cr significance from a UkeUhood-ratio 
test, based on the difference between x^ values of the best-fit 
NFW and SIS models: Ax^ = Xsis.min " XNFw.min = 72.6 for a 
1 degree-of-freedom. 

5. DISCUSSION AND CONCLUSIONS 

We have developed a method for improving the statisti- 
cal precision of cluster mass profiles, combining independent 
weak-lensing distortion, magnification, and strong-lensing 
measurements. This extends recent weak-lensing work by 
Umetsu et al. (201 1) to include the central strong-lensing in- 
formation in a stacking analysis, for full radial coverage. Our 
methods take into account the cosmic covariance from uncor- 
rected large scale structure projected along the line of sight 
(Hoekstra 2003; Hoekstra et al. 201 1), as well as the effect of 
different cluster redshifts, so that error propagation in terms 
of lensing efficiency of individual clusters can be properly av- 
eraged. 

We have applied our method to a sample of four sim- 
ilarly high-mass lensing clusters (A1689, A1703, A370, 
and C10024+17), for which we have previously identified 
multiply-lensed images and measured weak magnification 
and distortion effects from deep HST and Subaru observa- 
tions (Broadhurst et al. 2005a; Umetsu 8c Broadhurst 2008; 
Umetsu et al. 2010; Zitrin et al. 2010; Medezinski et al. 2010, 
201 1 ; Umetsu et al. 201 1). For our sample of massive clusters 
the strong and weak lensing regimes contribute equal loga- 
rithmic coverage of the radial profile and can be compared 
for consistency in the region of overlap. We have formed an 
averaged radial mass profile {T,{R)) from stacking the clus- 
ters (Figure 1), which shows a progressive steepening with 
increasing radius from R = 40kpc /;"' to 2800 kpc . The in- 
ner radial boundary is chosen to be sufficiently large to avoid 
smoothing from cluster miscentering effects (Johnston et al. 
2007), where the typical offset between the BCG and the dark 
matter center is estimated as d ^ 20 kpc /j"^ for our sample 
from our detailed strong-lens modeling (see Section 3.3). The 
stacked full mass profile is detected at a high significance level 
of 58c7 over the entire radial range. It is found here that ig- 
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Table 1 . Cluster sample and lensing data 



Cluster 


Redshift 


Einstein radius 


Strong lensing 


Weak lensing 


S/N 




z 




psl nsl 

min ' "max 




pwl nwl 
min' max 










(") 


(kpcr') 




(kpc/j"') 






A1689 


0.183 


53±3"(z. = 3.04) 


40, 125 


12 


129,2325 


11 


35 


A1703 


0.281 


31 ±3"(z. = 2.627) 


40, 177 


14 


179,2859 


10 


29 


A370 


0.375 


37±3"fc = 2) 


40, 149 


15 


152,3469 


14 


29 


C10024+17 


0.395 


30 ±3"(z. = 1.675) 


40, 126 


14 


134,3359 


12 


26 



Note. — For each cluster a joint mass profile is defined in N — N^^ discrete radial bins over the radial range 




^ 100 1000 

Projected radius, R [kpc//i] 

Fig. 1. — Top: the average projected mass profile S(R) (filled squares) with its statistical Icr uncertainty as a function of the projected radius R, which is 
obtained by stacking individual full mass profiles (thin gray lines) of four high-mass clusters (A1689, A1703, A370, and C10024-I-17 with My,,. > IO'^Mq at 
(z/) = 0.32) derived from Hubble strong lensing (R ^ 150 kpc /!"') and Subaru weak lensing (R ^ 150kpc /i"' ) measurements. The stacked mass profile exhibits 
clear continuous steepening over a wide range of radii, from R = 40 kpc to 2800 kpc /i"' 1.4rvii , which is well described by a single NFW profile (solid 
hne). The dashed line shows the contribution to the variance from uncoiTelated large scale structure projected along the line of sight. Bottom: the logarithmic 
slope of the stacked mass profile (open squares with error bars), d\n{T,) /dlnR, is shown as a function of projected radius along with the NFW model (solid line) 
shown in the top panel. The projected logaiithmic slope shows a clear continuous steepening with increasing radius, consistent with the NFW model. 



noring the cosmic noise contribution will underestimate the 
errors by ^ 30%-40%. This is due to the correlation of this 
noise between radial bins and can only be reduced by averag- 
ing over independent lines of sight, with uncorrelated line of 
sight structures, i.e. by averaging over well separated clusters. 

Our stacked projected mass profile with a continuously 
steepening radial trend is very accurately described by the 
NFW form predicted for the family of CDM-dominated ha- 
los, whereas it strongly disfavors the SIS model at 62a signif- 
icance. In the context of an assumed gNFW profile, the cen- 



tral cusp slope is constrained as a = 0.89^q 3^ (at r ^ 0.02rvii-; 
see Figures 1 and 2), being consistent with, but slightly shal- 
lower than, the simple NFW form with a= I. Our results are 
in agreement with recent high-resolution simulations, which 
find asymptotic inner slopes somewhat shallower than unity, 
TbdC?" ^ 0) ^ 0.9, for galaxy- and cluster-sized ACDM ha- 
los (e.g., Merritt et al. 2006; Graham et al. 2006; NavaiTO 
et al. 2010). Note NFW define this profile for halos which 
they identify as in virial equilibrium, in terms of the simulated 
CDM particles (see Section 2.2.2 of Navarro et al. 1997). The 
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Fig. 2. — Constraint on the gNFW model parameters, namely, the central cusp slope a, the halo virial mass Myi,, and the halo concentration f_2 = Cvir/(2 — «), 
when all of them are allowed to vary, derived from the averaged radial mass profile of A1689, A1703, A370, and C10024+17 shown in Figure 1. The left and 
right panels show the two-dimensional marginalized constraints on (Myi, , a) and (c_2, a), respectively. In each panel of the figure, the contours show the 68.3%, 
95.4%, and 99.7% confidence levels, and the cross indicates the best-fit model parameters. 



clusters we have selected for our stacked analysis are, in terms 
of their lensing properties, very well behaved with at most 
only ^ 10% perturbations in mass visible locally in the two- 
dimensional mass distribution, and otherwise very symmet- 
ric over most of the radius (Broadhurst et al. 2005a,b, 2008; 
Umetsu et al. 2010). Detailed hydrodynamical simulations 
show that equilibrium is relatively rapidly achieved in only a 
few sound crossing times after a major merger, though some 
dynamical and gas disruption may continue for over a Gyr 
This is not important in terms of the central relaxation time of 
the dark matter (Ricker & Sarazin 2001; Umetsu et al. 2010). 

An accurate measurement of the cluster mass profile en- 
ables us to constrain dark matter models. Recently Woo & 
Chiueh (2009) examined in detail an extremely light bosonic 
dark matter (ELBDM) model (m ~ 10^^ eV) as an alternative 
to CDM in the context of nonlinear cosmic structure forma- 
tion. ELBDM with a de-Broglie wavelength of astronomical 
length scales, if it exists, may well be in a ground-state Bose- 
Einstein condensate and hence well described by a coherent 
wave function, which may naturally account for the perceived 
lack of small galaxies relative to the ACDM model (Klypin 
et al. 1999; Peebles & Nusser 2010). Woo & Chiueh (2009) 
showed that, irrespective of whether halos form through ac- 
cretion or merger, ELBDM halos can form steepening density 
profiles of the form similar to the standard CDM, but with 
perhaps a steeper central cusp slope of 730 ~ 1 .4 and a shal- 
lower outer slope of 730 ~ 2.5. During a merger between 
condensates interesting large-scale interference occurs which 
will differ markedly from standard collisionless CDM, and it 
will be important to explore this class of dark matter further 
via more extensive and detailed simulations for testing against 
accurate lensing profiles of both relaxed and merging clusters. 

The mean concentration for the four massive lensing clus- 



ters considered here is found to be Cyir = 7.68^g4Q (at a 
mean virial mass Mvh- = 1 ■54;!;q [q x 1O'^M0 h~^), which is ap- 
parently higher than the standard ACDM predictions evalu- 
ated at the mean redshift (z/) = 0.32 of our sample: Cvi, = 
4.5;!;[ q (the errors quoted represent a Icr lognormal scatter 
of ^[logjgCvir] = 0.11) for relaxed clusters derived by Duffy 
et al. (2008) from A^-body simulations based on the WMAP 
5-year data and Cvk « 4.4 by Klypin et al. (2010) from the re- 
cent Bolshoi ACDM A^-body simulation. More recent results 
with greater mass resolution based on four large A^-body sim- 
ulations (Bolshoi, MultiDark, Millennium-I and II) exhibit a 
complex mass and redshift dependence of the median con- 
centration, namely a flattening and upturn of concentration at 
very high mass and redshift (Prada et al. 201 1). Accordingly, 
their concentrations derived for cluster-sized halos (i.e., rare 
objects corresponding to high- cr peaks in the primordial den- 
sity field) are substantially higher than previous results based 
on smaller simulations. Interestingly, they find a concen- 
tration of Cvir ~ 7 for their most-massive relaxed halos with 
Mvir ~ IO'^Mq/;-' at z = (Figure 15 of Prada et al. 2011). 
A comparison between our results and the ACDM predictions 
(Duffy et al. 2008; Klypin et al. 2010; Prada et al. 2011) is 
given in Figure 3. 

An accurate characterization of the observed sample is cru- 
cial for any cluster-based cosmological tests. In the extreme 
case, those clusters identified by the presence of a giant arc 
represent the most lensing-biased population. Calculations 
of the enhancement of the projected mass and hence boosted 
Einstein radii (say, 6'Ein > 20") find a statistical bias of ^ 34% 
derived from A^-body simulations of the ACDM model (Hen- 
nawi et al. 2007). Semi-analytical simulations incorporating 
idealized triaxial halos yield a ^ 50% bias correction (Oguri 
& Blandford 2009). Applying a conservative 50% bias cor- 
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Fig. 3. — Joint constraints on the mass and concentration parameters (A/vir,fvii) for a sample of four liigli-mass lensing clusters (A1689, A1703, A370, and 
C10024+17) derived from their stacked full mass profile (E(fi)) (Figure 1), compared to ACDM predictions (Duffy et al. 2008; Klypin et al. 2010; Prada et al. 
201 1) in the Cyir-Mvir plane. The cross shows the best-fit NFW parameters, and the contours show the 68%, 95%, and 99.7% confidence levels (Ax* = 2.3,6.17, 
and 11.8). The A'-body predictions of Duffy et al. (2008), Klypin et al. (2010), and Prada et al. (201 1) are shown as solid curves, with Icr lognormal scatter (taken 
from Duffy et al. 2008) indicated by the shaded area. Also shown are the levels of selection and projection bias for a strong-lensing cluster population derived 
from W-body (34%; dotted line) and semi-analytical (50%; dashed line) simulations, where the prediction by Duffy et al. (2008) is taken as the reference of the 
comparison. 



rection, we find a discrepancy of about 1 .8(7 with respect to 
the ACDM predictions by the Duffy et al. (2008) model for 
relaxed clusters (see Figure 3). If this large bias 50%) 
is coupled to a sizable intrinsic scatter in concentration, esti- 
mated for the full halo population to be (T[logioCvii-] =0.11- 
0.15, then our measurements can come into line with standard 
ACDM. 

The results presented here are very favorable in terms of 
the standard explanation for dark matter, as collisionless and 
non-relativistic, interacting only via gravity, with a very pre- 
cise match between our composite mean mass profile, and that 
of the general form of the mass profile advocated for massive 
halos in virial equilibrium. The relatively high concentration 
we obtain for the averaged profile is consistent with previ- 
ous lensing work which similarly detected a concentration 
excess in the lensing based measurements for many individ- 
ual relaxed strong-lensing clusters (e.g., Gavazzi et al. 2003; 
Kneib et al. 2003; Comerford & Natarajan 2007; Broadhurst 
et al. 2008; Oguri et al. 2009). This possibly interesting ten- 
sion between cluster lensing observations and ACDM mod- 
els can be more definitively addressed with fuU-lensing data 



for new cluster surveys, such as CLASH", LoCuSS, Subaru 
Hyper Suprime-Cam, and XMM-XXL (Pierre et al. 2010), to 
meaningfully examine the Cyir-Myir relation over a wider mass 
and redshift range when applied to sizable samples of relaxed 
clusters. It is highly desirable to cover the full profile by 
combining accurate weak and strong lensing measurements, 
requiring several sets of multiple images over a wide range 
of source redshift, to obtain a meaningful model-independent 
inner profile and to add weak lensing with sufficient color in- 
formation to exclude the otherwise sizable dilution effect on 
the weak lensing signal from foreground and cluster mem- 
bers. The CLASH survey is in particular designed to gen- 
erate such useful data free of systematics in both the weak 
and strong regime, with first results for the substantial smaller 
mass cluster A383 with M^i, = 5.37!°;™ x 10'''Mo /j"' (Zitrin 
et al. 201 Id) showing similar behavior (cyir = 8.77;|{J 42). 
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